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Using first principles calculations, we determined the epitaxial-strain dependence of the ground 
state of the 1:1 SrCrOs/SrTiCU superlattice. The superlattice layering leads to significant changes 
in the electronic states near the Fermi level, derived from Cr t 2 g orbitals. An insulating phase 
is found when the tensile strain is greater than 2.2% relative to unstrained cubic SrTiCU. The 
insulating character is shown to arise from Cr t^g orbital ordering, which is produced by an in-plane 
polar distortion that couples to the superlattice d-bands and is stabilized by epitaxial strain. This 
effect can be used to engineer the band structure near the Fermi level in transition metal oxide 
superlattices. 
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Metal-insulator transitions (MIT) occur in many tran¬ 
sition metal oxides as a function of composition, temper¬ 
ature, pressure and epitaxial strain |1 . At the atomic 
scale, the mechanism for the metal-insulator transition 
depends on the physics of the insulating state. For 
Mott insulators, the relevant parameters are the tran¬ 
sition metal (TM) d -level occupation, which can be 
changed through doping, and the bandwidth, which can 
be changed by varying the TM-oxygen-TM bond angles. 
For band insulators, a gap can be opened at the Fermi 
level by a structural distortion that lowers the crystal 
symmetry. A symmetry-breaking distortion that lifts 
the degeneracy of a partially occupied state will always 
lower the energy; this is called a Jahn-Teller distortion 
or a Peierls distortion depending on whether the degen¬ 
erate states are localized atomic orbitals or extended 
bands. Polar distortions have also been shown to couple 
to states near the Fermi level, proving useful in band-gap 
engineering |2j [3 . When the distortion can be controlled 
by an applied field or stress, the system can dynamically 
be driven back and forth through the MIT, a property of 
particular interest for applications to high-performance 
switching devices[4;, S\. 

SrCrOs (SCO) is a d 2 perovskite. It was ini¬ 
tially reported as a paramagnetic and metallic cubic 
perovskite |6 . Experimental studies also showed a tetrag¬ 
onal C type antiferromagnetic (C-AFM) phase, with 
the Neel temperature below 100K and space group 
P4/mmm, coexisted with the cubic phase in the low 
temperature [7, 8]. This tetragonal C-AFM state is fur¬ 
ther investigated in first-principles studies [9l [10] • It is 
correlated to a partial orbital ordering 
where the d xy is nearly occupied while d yz and d xz or¬ 
bitals are both half occupied j8j 9.. For the single layer 
film of SCO epitaxially grown on SrTiOs, due to the miss¬ 
ing apical oxygen of the CrOs octahedra at the surface, 
d yz and d xz orbitals become lower in energy than d xy , 
leading to a d yz dl, z d^ y orbital ordering and an insulating 
state [TT|. However, recently Zhou et al claimed that SCO 


is an insulator while the insulator-metal transition occurs 
under sufficient pressure [12], due to the bond instability 
found around 4 GPa, suggesting the strong interaction 
between structure and electronic structure. 

Layering a d-band perovskite oxide in a superlattice 
with a d° insulating perovskite will lead to large changes 
in the d-band states near the Fermi level, depending on 
the thickness of the d-band layer fl3fil5] . For SCO, 
SrTiOs (STO) is a natural choice as the second com¬ 
ponent. It is a d° insulator with a 3.25 eV band gap and 
an A-site cation in common with SCO. In addition, the 
structure of STO can be tuned by epitaxial strain, with 
different octahedral rotation patterns and polar distor¬ 
tions in compressive and tensile strain pRHl8] : these dis¬ 
tortions are expected to influence the structure in the 
SCO layer. 

In this Letter, we present first-principles results for the 
ground-state phase sequence of the 1:1 SrCrOs/SrTiOs 
superlattice with varying epitaxial strain, with a first- 
order MIT from a metallic x-type P2\je phase to a polar 
insulating G-type Pmm2 phase observed at 2.2% tensile 
epitaxial strain. We show that the insulating character 
arises from orbital ordering induced by a polar distortion 
in the SCO layer. The polar distortion is shown to be 
stabilized by a combination of epitaxial strain and the 
polar distortion in the adjacent STO layer. This offers 
the possibility of driving the system through the MIT by 
an applied field or stress that coupled to the polar mode. 

Our first-principles calculations were performed us¬ 
ing the local density approximationpLQj [20] with Hub¬ 
bard U (LDA+U) method implemented in the Vienna 
Ah initio Simulation Package (VASP-5.2j21, 22]). We 
used the Dudarev implementation 123^ with effective on¬ 
site Coulomb interaction U — 1.5 eV to describe the lo¬ 
calized 3d electron states of Cr atoms. The projector aug¬ 
mented wave (PAW) potentials [24j [25] contain 10 valence 
electrons for Sr (4s 2 4p 6 5,s 2 ), 12 for Cr (3p 6 3d 4 4s 2 )10 for 
Ti (3p 6 3d 2 4s 2 ) and 6 for O (2s 2 2p 4 ). x a/2 x 2 and 
2x2x2 cells were used to allow the structural dis- 
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tortions and magnetic orderings for the 1:1 superlattice. 
500 eV energy cutoff, 4 x 4 x 4 k mesh and 5 x 10 _3 
eV/A force threshold were used for structural relaxations. 
A 15 x 15 x 11 k-point mesh was used for density of 
state calculations. Our calculations give a = 3.812A for 
the relaxed tetragonal state of C-AFM SrCrOs, consis¬ 
tent with the previous first-principles and experimen¬ 
tal studies[H [9], with aspect ratio c/a = 0.964 slightly 
smaller than the experimental value of 0.99. The epi¬ 
taxial strain value is defined with respect to 3.849A, the 
computed lattice constant of the cubic STO. The effects 
of epitaxial strain were studied through “strained bulk” 
calculations [26]. 

At each value of epitaxial strain in the range —4% - 
3%, we determined the ground state (GS) structure of 
the superlattice using the “stacking method” [18,. This 
involved combining the computed low-energy structures 
of epitaxially strained STO and SCO to obtain a set 
of starting structures, which were then relaxed to the 
nearest energy minimum using first-principles calcula¬ 
tions. To obtain the low-energy structures of SCO, we 
considered F, G, A, and C-type AFM ordering in the 
tetragonal phase with various distortions. C-AFM was 
found to be the ground state magnetic ordering in the 
full range of epitaxial strain considered. For C-AFM 
and FM states, PA/mmm is the GS structure. For G- 
AFM, PA/mmm is the GS from —4% to 0%, P4/mbm 
Jahn-Teller (JT) [27l [28] distortion is the GS beyond 0% 
and anti-polar distortion Pmcm is also unstable beyond 
2%. For A-AFM, PA/mmm is the GS from —4% to 
0%, PA/mbm Jahn-Teller (JT)[27, 28> distorted struc¬ 
ture is the GS beyond 0%, and polar distortion Pmm2 
is the unstable beyond 3%. The low-energy structures of 
epitaxially-strained STO have been previously discussed 
in the literature [T6HT8] . We use the low-energy struc¬ 
tures of STO listed in Table I of Ref. [18 i , consistent 
with the reports in other studies. When the epitaxial 
strain is within —1% - 1%, the low-energy structures of 
STO are generated by the M and R point octahedral rota¬ 
tions. For compressive epitaxial strain greater than —1%, 
the out-of-plane polar mode becomes unstable, while the 
in-plane rotations are suppressed. For tensile epitaxial 
strain greater than 1%, the main distortions are in-plane 
polar modes along [110] and [100] and out-of-plane rota¬ 
tions are suppressed. 

To include magnetic ordering in the stacking method 
for the 1:1 superlattice, we consider three magnetic 
phases for the single layer of SCO that can fit in the 
2x2 lattice, as shown in Fig. [lja). The three phases are 
FM, G-AFM for the checkerboard ordering in the CrO 
plane, and x-type AFM for the AF ordering along [100] 
and FM along [010] in the CrO plane. For each magnetic 
state of SCO and each value of epitaxial strain, we gen¬ 
erate the starting structures by following the procedures 
of the stacking method. We also carry out relaxations 
for a small number of random starting structures p~8] . 
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FIG. 1. (a) Magnetic orderings in the SCO layer. Red: fer¬ 
romagnetic (F), black: G-AFM (G), and blue: x-type AFM 
(x) states, (b) GS structures and energies for F (red), G 
(black), and x (blue) magnetic states as functions of epitaxial 
strain. The solid curves guide the eye. Insulating and metal¬ 
lic states are denoted by open and solid symbols. Shapes of 
data points indicate the space groups: pentagons ( Pbam ), 
triangles (P2/ra), circles (P4/ra&ra), diamonds (P2i/c), five- 
pointed stars (Pc) and squares (Pmm2). The inset shows the 
energies of low energy structures relevant to the F state in the 
epitaxial strain range 1.8% to 2.9%. 

Thes show that from —3% to —1% epitaxial strain, the 
JT distortion is found in the SCO layer, while it is not 
metastable in the bulk state. 

In FigjTJb) we plot the total energies of the GS of the 
superlattice obtained by the “stacking method” for the 
three magnetic states versus epitaxial strain. From —3% 
to 2.2%, the x-AFM is the GS magnetic state, as expected 
from the resemblance of magnetic ordering to the bulk C- 
type AFM state. Beyond 2.2%, the GS magnetic state 
changes to G-AFM. 

It is instructive to examine the epitaxial strain depen¬ 
dence of the structure for each magnetic ordering. For 
FM state, from —3% to 0%, PA/mbm structure with the 
a°a°c~ oxygen octahedral rotation pattern is the lowest 
energy structure in this phase. From 0% to 2.3%, the 
lowest energy structure is P2i/c, with an a~a~c~ rota¬ 
tion pattern and JT Q 2 mode. The amplitude of the Q 2 
mode is about 1% of the lattice constant and does not 
change much with the strain. As the tensile epitaxial 
strain increases, the a°a°c~ rotation decreases while the 
a~a~c° increases. Beyond 1%, the lowest energy struc- 













3 



FIG. 2. The PDOS of the Cr d yz orbital for G-type 1:1 
SCO/STO superlattice and C-type bulk SCO. 


ture of FM state is insulating. Beyond 2.3%, the lowest 
energy structure is Pc, with the a~a~c° rotation and 
small and decreasing amount of a°a°c~ rotation, JT Q 2 
mode in the SCO layer, and polar distortions along [110]. 

For the G-AFM state, from —3% to —1%, Pbam is the 
lowest energy structure, with the a°a°c~ rotation mainly 
in the STO and JT Q 2 distortion in the SCO layer. The 
JT distortion amplitude is small but enough to lift the 
d yz d xz degeneracy. However, the system is still metallic 
due to the band overlap. From 0% to 2.2%, P2\/c with 
the a~a~c~ oxygen octahedral rotation pattern and JT 
Q 2 mode is the lowest energy structure. Beyond 2.2%, 
the polar structure Pmm2, with only polar distortion 
along [100] becomes the lowest energy structure, leading 
to the insulating overall GS, with magnetic moment 1.8 
Pb- 

For the x-AFM state, from —3% to 0%, the lowest en¬ 
ergy structure is nonpolar P2/m, with a°a°c~ rotation 
mainly in the STO layer, and small JT Q 2: Q 3 modes in 
the SCO layer. From 0% to 2.2%, the lowest energy struc¬ 
ture is P2i/c. Beyond 2.2%, the lowest energy structure 
is insulating. In particular, from 2.2% to 2.7%, the low¬ 
est energy structure is Pc with the a~a~c° rotation and 
a small and decreasing amount of a°a°c~ rotation, JT 
Q 2 mode in the SCO layer and polar distortion along 
[110]. The in-plane rotation is now the main distortion 
pattern, and the amplitudes for both STO and SCO lay¬ 
ers are similar, due to the connection of oxygen octahe- 
dra. Beyond 2.8%, the lowest energy structure changes 
to Pmm2. 

The phase boundary at 2.2% epitaxial strain is of par¬ 
ticular interest. The inset of Fig. [l] shows the energies 
for F, G, and x structures for 1.8% — 2.9% epitaxial 
strain. The energy of the G-type polar Pmm2 struc¬ 
ture decreases relative to that of the x-type nonpolar 
P2\/c structure with increasing epitaxial strain, with a 
first-order transition at 2.2% epitaxial strain. We note 
the metal-insulator transition associated with the change 
from the nonpolar to the polar structure. 



FIG. 3. PDOS of the spin up Cr t 2g in the 1:1 superlattice for 
+3% strain, with G-AFM magnetic ordering. The left panels 
represent the PDOS without (top) and with (bottom) the 
polar distortion. The vertical dashed line marks the energy 
of the highest occupied state. The distortions in the SCO 
layer are shown on the right part, where gray, blue and orange 
spheres represent Sr, Cr and O ions, respectively. 


To figure out the mechanism of the insulating phase 
in the 1:1 superlattice beyond 2.2% epitaxial strain, we 
first consider the effect of the superlattice layering on the 
bands near the Fermi level. It eliminates the dispersions 
of Cr t 2g bands along k z , and significantly narrows the 
band widths of d yz and d xz . Fig. [2] shows the PDOS of 
d yz of one Cr atom in the G-type 1:1 superlattice and 
C-AFM bulk SCO for 3% tensile epitaxial strain in the 
PA/mmm symmetry. The 1:1 superlattice is metallic for 
the PA/mmm structure, with the d yz band much nar¬ 
rower than in the bulk case. The thickness of the STO 
layer has little effect. Our calculations show that the Cr- 
layer structure and bands are very similar for 1:1 and 1:3 
superlattices. 

Next, we consider the coupling effect of the polar dis¬ 
tortion to the superlattice insulating Pmm2 state. The 
Pmm2 structure is generated from the high symmetry 
reference structure by the doubly-degenerate in-plane 
E u polar distortion, with eigenvector (0.00, 0.00, —0.15, 
0.08, 0.53, 0.53, 0.37, 0.40, 0.06, 0.33) specifying the dis¬ 
placement pattern of the atoms (Sr, Sr, Cr, Ti, O^i, O z2 , 
0,1, 0,2, O y i, O^)- From 2.2% to 3%, the band gap 
increases from 0.27 to 0.47 eV, and the polarization in¬ 
creases from 36 to 41 fiC/cm~ 2 . 

The effects of the in-plane polar distortion on the elec¬ 
tronic structure of the superlattice are evident from the 
projected density of states (PDOS) of the t 2g bands of 
the spin up Cr atom, shown in Fig. [3] The layered struc- 
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FIG. 4. The evolution of PDOS of Cr d xy , d yz and d xz as a 
function of polar displacement of Cr ion in [100]. The ampli¬ 
tude of the displacement for each curve is given in the legend 
in units of the in-plane lattice constant. 

ture of the superlattice splits d xy from d xz and d yz and 
increasing tensile epitaxial strain lowers the energy of the 
d xy orbital relative to d yz and d xz . At 3% epitaxial strain, 
in the undistorted structure, all three d orbitals are par¬ 
tially occupied. The in-plane polar distortion lifts the 
degeneracy of d yz and d xz , so that in the polar Pmm2 
state the d xy and d yz orbitals are fully occupied while d xz 
is unoccupied, corresponding to d] cy d yz d^ z orbital order¬ 
ing. 

In fact, the 00 can be produced by polar displace¬ 
ments just of the Cr atoms. To show this, we consider a 
Pmm2 1:1 SCO/STO superlattice, in which only the Cr 
atoms are uniformly displaced along [100], while all other 
atoms stay at the high symmetry positions. In Fig. [4]we 
show the PDOS of the distorted superlattice as a func¬ 
tion of Cr displacement. There is a dramatic downward 
shift of the d yz with increasing displacement, accompa¬ 
nied by a smaller upward shift of d xy and a decrease in 
the d xy bandwidth. 

Finally, we discuss why the in-plane polar distortion 
becomes the GS for tensile strain. Our calculations show 
that for any distortion that lifts the d yz and d xz degen¬ 
eracy, there is a linear energy gain. This explains the 
emergence of JT distortions in P2\jc or Pbam struc¬ 
tures for the 1:1 superlattice. However, in general, for 
large tensile epitaxial strain, the in-plane polar state is 
likely to be more favorable than the JT distortion due to 
the well-known polarization-strain coupling. The polar 
distortion in the STO layer is also an important influ¬ 
ence. To show this, we fixed the STO layer in a nonpolar 
structure for 3% strain, froze in the in-plane polar mode 


and JT distortion in the Cr02 layer in turn, and relaxed. 
Distortions survive in both cases and the relaxed struc¬ 
ture with JT distortion has lower energy, suggesting that 
the in-plane polar mode in the STO coupled differently 
to the in-plane polar mode and the JT mode in the SCO 
layer. 

This insulating state also raises the possibility of con¬ 
trolling band gap by applied electric field. Given the 
SCO/STO superlattice in the insulating state, an in¬ 
plane electric field will change the atomic positions, and 
hence change the band gap, because the band gap is sen¬ 
sitive to the displacement of the Cr atom relative to the 
O atoms around it. 

In summary, we have studied both the lattice and 
electronic structures of the ground state for the 1:1 
SCO/STO superlattice. Distortions in SCO layers are 
established by the superlattice layering with STO. For 
tensile epitaxial strain, due to the in-plane polar dis¬ 
tortion associated with nonzero Cr displacements, the 
d xy dy Z d xz orbital ordering can be formed and the band 
gap is therefore induced. The polar distortion induced 
metal-insulator transition can be used to engineer the 
SCO band structures near the Fermi level. Our study 
sheds light on a new way to control electronic band struc¬ 
tures and approach the metal-insulator transition point. 
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